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A variety of observational tests seem to suggest that the universe is anisotropic. This is incom- 
patible with the standard dogma based on adiabatic, rotationally invariant perturbations. We point 
out that this is a consequence of the standard decomposition of the stress-energy tensor for the 
cosmological fluids, and that rotational invariance need not be assumed, if there is elastic rigidity 
in the dark energy. The dark energy required to achieve this might be provided by point symmetric 
domain wall network with P/p — —2/3, although the concept is more general. We illustrate this 
with reference to a model with cubic symmetry and discuss various aspects of the model. 
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PACS numbers: 

Observations of the angular power spectrum of the 
cosmic microwave background (CMB) at intermediate 
scales, multipoles 50 < £ < 1000, made by the WMAP 
satellite and other experiments combined with obser- 
vations of type la supernovae at high rcdshift |^ Q 
and the galaxy power spectrum 0,0 appear to have con- 
firmed the basic tenets ACDM model based on adiabatic 
initial conditions created during inflation. However, ob- 
servations on larger scales, while nominally compatible 
with scale-invariance, and hence with the ACDM model, 
appear to be incompatible with the assumptions of Gaus- 
sianity, isotropy, or both 

Of particular interest here are the North-South asym- 
metry in estimates of the power spectrum in the direc- 
tion of {I, b) « (57, 10) and the so-called "axis-of-evil" in 
the derived multipolc coefficients in a coordinate system 
orientated such that the North pole is in the direction 
of {l,b) « (—110,60), which is orthogonal to the North- 
South asymmetry. These results have led T. Jaffe et al 13 
to consider the possibility that the observed CMB has 
two sources: a Gaussian component based on approxi- 
mately scale invariant adiabatic fluctuations, and a tem- 
plate created from a low density Bianchi Vllh universe. 
They found that such a model with the template oriented 
with a rotation axis in the direction {l,h) ~ (222, —62), 
shear cf/Hq = 2.4 x 10^^° and right-handed vorticity 
uj/Hq = 6.1 X lO^^*' fitted the observed anisotropics. 
This result is independent of frequency, which appears 
to rule out an obvious galactic origin for this effect and, 
therefore, suggests an anisotropic universe which is also 
rotating. There is also some evidence for these phenom- 
ena in the COBE data but with lower signal-to-noise. 

While this work illustrates a fundamental aspect of the 
data, the two components used in this analysis are log- 
ically incompatible from theoretical point of view: the 
first is based on a flat universe created during inflation 
and the second is an anisotropic universe with sub-critical 
matter density. In more recent work 0, it has been 
shown that one cannot create a sufficiently strong effect 
to explain the data when the universe is flat and domi- 
nated by a cosmological constant. This begs the question 
which we attempt to address in this Letter : can there 
be a cosmological origin of the observed anisotropy and 



can it be achieved within a dark energy/CDM model? 
A possible solution, we suggest, lies in the realization 
that the stress-energy tensor of the cosmic fluids need 
not be rotationally invariant at pcrturbative order, and 
that the standard scalar-vector-tcnsor (SVT) split of the 
linearized gravity and conservation equations can only 
be made when the full set of isometrics of 3D Euclidean 
space are assumed. This need not be the case if the dark 
energy component, rather than being a cosmological con- 
stant, is described by th e dy namics similar to those of an 
elastic continuum solid [lOl llll | which might have a mi- 
croscopic realization in a static domain wall conflguration 
where the density, p, and the pressure, P, are related by 
P/p ~ —2/3 and there is sufficient rigidity 0, 0| to 
achieve stability. The current constraints on such a sce- 
nario are presented in ref. [T^ under the assumption of 
isotropy. 

The standard lore of cosmological perturbation theory 
(see, for example, ref. |l5j ) is to decompose the perturbed 
stress-energy tensor for each of the cosmological fluids, 
5T^^, into eigenfunctions of the rotationally invariant 
Laplacian. In particular, one typically writes 5T^i, — 
{Sp + 5P) u^^u^ - (5P(5'^^ -f (p P){y>'Uy -h u^'Vy) -h n''^, 
where is a unit vector specifying the fluid flow lines 
{u^Uv = 1), (5/9 is the density perturbation, 6P is the 
pressure perturbation, (u^V^ = 0) is the velocity per- 
turbation orthogonal to the flow and 11'^,^ (II'^^u^ = 0, 
W^ijvy ~ and 11^^ = 0) is the anisotropic stress per- 
turbation. If one transforms the spatial coordinates to 
Fourier space with wave- vector ki = kki, then one can 
perform the SVT split for the velocity perturbation as 
Vi = V^ki+V^^li+V^'^rhi and the anisotropic stress per- 
turbation as n^j = {k^kj-Sij /3)n^+ki{U^^ij+U^'^Thj)+ 
kj (n^i -t- U^'^ihi ) + n+ {iilj -■m,mj) + W {krhj + rhjj ) , 
where U and rhi are unit vectors which form an orthonor- 
mal triad with fcj. is the irrotational velocity and 
V^^ and V^"^ are the components of the vorticity. 11^ is 
the scalar anisotropic stress, 11^^ and 11^^ are the vec- 
tor anisotropic stresses and 11+ and 11^ are the tensor 
anisotropic stresses. 

The metric perturbations can also be split in a similar 
way and due to rotational invariance the perturbed equa- 
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tions of motion can be split up into three non-interacting 
blocks, the scalars, vectors and tensors, which can be 
evolved separately. In the case where the initial con- 
ditions are density waves (scalars), then the vorticity 
(vectors) and the gravitational waves (tensors) will re- 
main zero, and the statistical rotational invariance and 
the Gaussianity of the initial fluctuations will be main- 
tained at linear order. We stress that this split, while 
well motivated, particularly in the context of describing 
the radiation and cold dark matter (CDM) components, 
is only an assumption. We will show that in a more 
general treatment of cosmological fluids, as might be re- 
quired to describe dark energy, one need not make this 
assumption. 

Except for the case of a cosmological constant, it is 
necessary to consider the evolution of perturbations in 
all other fluid based dark energy models since it is im- 
perative to maintain energy conservation. This was first 
considered in the context of scalar field dark energy in 
ref. and a number of recent works (for example, 
refs.|l7l lisll') have developed the treatment of this issue. 
Scalar field fluids are not adiabatic since the perturba- 
tions in the scalar field allow for the CDM and dark en- 
ergy to have different rest frames. If one wishes to con- 
sider general adiabatic fluids, then one has to consider 
fluids whose macroscopic Lagrangian, >C, is just a func- 
tion of the metric. In this case the stress-energy tensor 
is given by the appropriate functional derivative of the 
Lagrangian with respect to the metric and one can also 
define a rank four tensor, W^^'^f"^ , which is the second 
functional derivative of the Lagrangian, such that 

T^'^ = -2\g\-'^'^^{\g\^'^C), (1) 
= ^2\g\-^/^-^{\g\^/^Tn, (2) 

dgpa 

where \g\ is the determinant of the metric. 

One can rewrite (0) as ST^"^ = 
-^{W^P" +Tf"'gP'')6gpa, where the perturbation 
in the metric is given in terms of the background metric 
perturbation hp^- and the Lagrangian perturbation 
of the fluid by Sgp^ = hp^ + 2V(p^cr)- Since the 
flow is defined relative to constant density lines, then 
one can decompose ^3 the stress energy tensor as 
TP'^ = puPu'^ + P'"' and the second functional derivative 
as W^'-'^P" = EP'^P" + pP-^uPu" + PP'^uP-u" - PP-Pu^u" - 
pP^u^uP - P^PuPu" - P'"^uPuP - puPu^uPu", where PP" 
[PP'^u^ — 0) is the pressure tensor and EP'^p'^ (which 
satisfies EP'^p'^Uc = and EP'^p'' = E^p^^^p"^^ = Ep^p") 
can be interpreted as an elasticity tensor |23| which 
in general has 21 components. One of these, the bulk 
modulus, is specified by the pressure and the other 20 
are shear moduli. It is this tensor which will give us 
a general parameterization of linearized perturbations 
in these adiabatic dark energy models. A detailed 



exposition of the isotropic case, where there is just a 
single shear modulus, is presented in ref. [U, suffice 
to say that stability requires the shear modulus to be 
sufficiently large to overcome the natural instability of 
fluids with negative pressure. 

Here, we will be concerned with anisotropic models 
where, by analogy to the standard theory of elasticity 
in solids, we can deduce that all the possible cases can 
be completely classified in terms of the well-known Bra- 
vais lattices. In particular, the perturbations must have 
cubic, hexagonal, rhombohedral, tetragonal, orthorhmo- 
bic, monocline, or triclinic symmetry |23|. In the rest 
of this work we will consider for definiteness the cubic 
case for which the pressure is isotropic PP'^ = PjP'^, 
where 7^1/ — g^u — u^u^, and there are two non-zero 
shear moduli ^l, fJ-T plus the bulk modulus defined by 
/3 = (p + P)^- If one defines 1 = xx, 2 = yy, 3 = zz, 
4 = yz, 5 = xz, 6 = xy then the non-zero components 
of the elasticity tensor are given by E^^ = E'^'^ = E^^ = 
(3 + P + Ipl, £^12 ^ E^^ = E^^ = 13 - P - and 
^ ^55 ^ ^ p jf ^ ti-^ig returns to 

the isotropic case considered in ref. |21|. Moreover, moti- 
vated by domain walls, we will concentrate our numerical 
work on the case oi w = P/p = —2/3. The basic qual- 
itative features of our analysis will be present in more 
general cases. 

From the point of view of the present discussion, 
the important aspect of these cubic models is that the 
three independent speeds of wave propagation, com- 
monly called sound speeds, depend of the direction of 
the wave-vector, defined here in terms of the polar an- 
gles 6 and 0, as well as its amplitude ^3 as is the 
case in crystals. By requiring that each of the sound 
speeds are greater than zero for all directions, we showed 
that the lattice will be stable to continuum modes if 
both fih/p > 1/6 and pt/p > 1/6. Moreover, assum- 
ing Nambu-Goto walls, we have computed the two mod- 
uli for the three primitive cells with cubic symmetry, 
the Wigner-Seitz cells of simple cubic, body-centre cubic 
(BCC) and face-centred cubic (FCC) which correspond 
to polyhedral cells made from cubes, truncated octahc- 
dra and rhombic dodecahedra. The stability conditions 
are violated by the BCC cell, and the other two lattices 
have zero modes, meaning that the sound speed is zero in 
at least one direction. In this work we will assume that a 
stable structure can be constructed from a compound (as 
opposed to primitive) cell and discuss the phenomenology 
of the resulting perturbation equations for stable lattices. 

The equation of motion for (using a gauge where 
£,Pup = 0) is modified from that presented in ref. to 

[p + P)(f + HC) - 3f3He - Pid'dj^' + d'h/2) 
-p^id^djC + + d^h'j - d'h/3) ^ ApF\ (3) 

where Ti. is the conformal time Hubble parameter and we 
have used the synchronous gauge, h is the trace of the 
spatial metric perturbation hij , and 
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\ (9,9^ + dydy)e + + dy^y) + d^h\ + dyh% J 



The degree of anisotropy is quantified by A/i = — fJ-h 
which is zero in an isotropic model. If both /il and /ix 
are zero then this equation describes perturbations in a 
perfect fluid. 

If = 0, one can perform the standard SVT split 
since the isometrics of Euclidean space allow one to define 
pure SVT modes. However, if one attempts to perform 
the same split in the more general case, then the SVT sec- 
tors, which are usually decoupled, can source each other 
and initial conditions which comprize of pure scalar, adi- 
abatic modes can excite vector and tensor modes sponta- 
neously. In order to illustrate this wc have computed the 
power series solution to the equations of motion ^ cou- 
pled to the relevant linearized Einstein equations, plus 
those for the perturbations to radiation, CDM and dark 
energy (with w = —2/3) components, with initial condi- 
tions = Cde = and h,j = &k-^/'^{k,kj - \5rj). The 
first term in the expansion for is 




Bkx ky 
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Bky hz 
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(5) 



where A = - /ii), B = -A/i/2, fii^ = Ml/p, At = 
/ix/p and A/t is similarly defined. 

One can define the "would be" scalar displacement 



^DE 



^i^DE which is given to the same order by 
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(l-4AL)-12AA(fcp^ 



Hyl^Z 



(6) 

The first term, which is independent of the direction of ki , 
is what was computed in the isotropic case , but the 
second term, which is explicitly symmetric under cubic 
transformations, is direction dependent. Moreover, the 
equivalent vector displacement is non-zero if A^ ^ 0. 

One can define the density contrast in the dark energy 
component to be Sue = —(1 + w)(fc^QE + ^h) and its ve- 
locity perturbation to be V^^ 



V2- 
DE ^ 



Ci _ T/S Li 
SDE — ''de'^ 



This allows the definition of the total density 
contrast (5x = r^DE&E + f^miJm + ^^ri^r in terms of the den- 
sities of the dark energy, CDM and radiation components 
relative to critical, which is responsible for the gravita- 
tional potential and hence all observational effects, be 
they in the cosmic microwave background (CMB) or the 
galaxy distribution. From this we can then define the 
total power spectrum Px = |(5xP- 

We have evolved the equations of motion numeri- 
cally for w = -2/3, fLL = 0.18 and A/i = 0.01. 
Fig. n shows that time evolution of the scalar and vec- 
tor velocities of the dark energy component, |V[fg| and 





FIG. 1: Perturbation evolution for k = 10~^Mpc~^ for a 
component with w = —2/3, fiL = 0.18 and Afi = 0.01 in the 
direction 6 = Tr/2 and </> = n/8. On the left is the scalar 
velocity of the dark energy component |VdeI (solid line) and 
the vector velocity component (|VdeP + I^ideP)^''^ and on 
the right are the vector (dotted line) and tensor metric com- 
ponents. Note that the vector and tensor perturbation are 
non-zero even though the initial conditions were pure scalar. 
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FIG. 2: Pt (left) and \Vii\'' + \V^i\^ (right) at fc = 
lO^'^Mpc^^ as function with 6 and 4> plotted in the Hammer- 
Aitoff projection when w — —2/3, /11 ~ 0.18 and A/i — 0.01. 
Note the anisotropy of the Pt and that (I'Fde |^ -I- IV^ip)^''^ 
is non-zero. Both have obvious cubic symmetry. 



(I^deI^ + I^DEP)^^^ foi- k ^ lO-^Mpc"^ in the direc- 
tion 6* = 7r/2 and (p = tt/8, and also the vector and 
tensor metric components. The important point is that 
the vector velocity and these metric perturbations would 
be zero if A/t = 0, but are clearly non-zero here illustrat- 
ing the mixing of SVT modes. In Fig. |21we present the 



angular distribution of Px and IV^qeI"^ + I^4)eI^ 
present day for k ~ 10^'^Mpc^^. It is clear that there 
are cubic anisotropics at the level of a few percent in the 
distribution of Px, and that the vector velocity is highly 
anisotropic, being zero in some directions and large in 
others. 

In order to quantify the amplitude of the effect as a 
function of A/i we have computed two average quantities 
as function of fc. First, the normalized variance of the 
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FIG. 3: The normalized variance of the power spectrum K 
(left) and the ratio of the vector and scalar velocities R 
(right) for w = —2/3 and jiL = 0.18. Curves show vary- 
ing levels of the anisotropy A/i — 10~^(solid), 10~^(dotted), 
10"^ (short - dash), 10~^(long - dash). 

power spectrum over the angular directions, which 
is related to the kurtosis of the density field. This de- 
termines the level of anisotropy, or non-Gaussianity 
expected. The second is the ratio of the vector and scalar 
velocities of the dark energy, i?, which quantifies the level 
of local rotation. Formula for K and R are given by 

(7) 

where (..) corresponds to the average over the sphere. 

These quantities are plotted in Fig. |21for range of val- 
ues of the Afi. We see that the level of anisotropy peaks 
at around k 10~^Mpc~^ and increases with Afi as one 
would expect. This is because, for the specific choice 
of /11 which we have used, the Jeans length of the dark 



energy component is a substantial fraction of the hori- 
zon. Therefore, when a particular mode comes inside the 
horizon it only grows for a short period of time and the 
anisotropy is maximal just inside the the present horizon, 
that is, on large scales and for low multipoles of the CMB. 
The ratio of the vector and scalar velocities is constant 
on very large scales, the amplitude again increasing with 
Afi, and falls off on smaller scales since the vector (vor- 
ticial) velocity component decays once the mode comes 
inside the horizon. 

Therefore, we have shown that in principle the specific 
model which arc discussing gives rise to the anisotropy on 
the very largest scales, primarily since the dark energy is 
dominating when these scales cross the horizon. Qualita- 
tively, this kind of phenomenon has been observed in the 
WMAP data. It is clear from the preceding discussion 
that the power spectrum measured in small regions will 
have an excess variance over an isotropic, Gaussian case. 
Moreover, due to the point symmetry, it is inevitable 
that different Fourier modes will be coupled together on 
large scales inducing something qualitatively similar to 
an "axis-of-evil" . We should emphasize that we do not 
claim at this stage that we have shown any quantitative 
agreement between the predictions of this model and the 
observed anomalies. This is the next step in our work 
where we will compute the correlation matrix (airnO-i>m')- 
Finally, we make the disclaimer that we have also not yet 
shown that such a domain wall lattice can be formed in 
any reasonable scenario although we believe it to be pos- 
sible. Nonetheless, we feel that the work presented here 
illustrated an interesting avenue for future investigation. 
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